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Necessary conditions for separability are most easily expressed in the computational basis, while 
sufficient conditions are most conveniently expressed in the spin basis. We use the Hadamard matrix 
to define the relationship between these two bases and to emphasize its interpretation as a Fourier 
transform. We then prove a general sufficient condition for complete separability in terms of the spin 
coefficients and give necessary and sufficient conditions for the complete separability of a class of 
generalized Werner densities. As a further application of the theory, we give necessary and sufficient 
^\ • conditions for full separability for a particular set of n-qubit states whose densities all satisfy the 

' Peres condition. 
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^ | The study of non-classical correlations has led to a number of suprising results arising from the existence of entangled 
. states of separated subsystems . This has led to renewed interest in the study of entanglement itself (||| as wen 
as in applications such as quantum information theory and quantum communication Q . Before one can effectively use 
entanglement, it is necessary to determine if a given state p actually has entangled subsystems. It is this "separability" 
problem with which we concern ourselves in this paper. 

There exists a useful, general necessary condition for separability || and a theoretical necessary and sufficient 
condition but no operational necessary and sufficient conditions, and as a result attention has tended to focus 
\q \ on classes of densities ]To|-[l^]. In this paper, we record a useful variant of the Peres (necessary) condition and prove 
y—{ . a new sufficient condition for full separability of mixed states of a system composed of n-qubits. To do that, we 
r—i ' highlight the roles of the computational basis, composed of projections and raising and lowering operators, and the 
spin basis, composed of the identity and the real Pauli matrices. We derive a change of basis formula which facilitates 
changing from one basis to another and apply these insights to obtain the general sufficient condition for separability 
and to obtain both necessary and sufficient conditions for a particular class of states satisfying the Peres condition. 
In a separate paper we will show how these ideas generalize to higher dimensional states. 

A state defined on the Hilbert space Ha! ® Ha 2 is said to be separable if it can be written as 



o\ 

S3 
a: 

p = p(a)p (a) = y p(a)p A ^ (a) ® (a), (1) 

a , 

where p (a) = p Al (a) ® p A2 (a) with p Ak (a) a state on Ti.A k , and the p(a) are positive numbers that sum to one. Peres 
showed that a necesary condition for a density matrix of such a bipartite system to be separable is for its partial 
transpose to be a density matrix ||, where the partial transpose p Tl of p is defined by (a b\p Tl \a' b') = (a' b\p\a b'). 
(p T2 is defined analogously.) For 2 ® 2 and 2 ® 3 systems this condition is also sufficient ||. The Peres condition is 
basis independent, but to facilitate applications we derive a weaker version which is most usefully expressed in the 
computational basis. This result is based on the positivity of the subsystem states, the assumed degree of separability 
and the Cauchy-Schwarz inequality. First, we introduce some notation. Let j denote an n-long binary index vector 
and let and 1 stand for the n-bit numbers consisting of all 0's and all X's, respectively. The binary complement of j 
will be denoted by j = l®j where the addition is mod2. We shall write j — j 1 ] 2 to mean that j is the concatenation 
of j 1 and j 2 . Now assume p has the form given in (|l|). Then 



VWp\f)VW\k) = 



J2p(a)Pj,j (a) 



1/2 



y,p{a) Pk,k (a) 



1/2 



P(a) \J pfiji i a )pt\k^ ( a ) v p P,P ^Pklk 2 ( a ) ( 2 ) 
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where the first inequality is the Cauchy-Schwarz inequality, the middle equality reflects the assumed separability, 
and the second inequality follows from the positivity of the subsystem density matrices. Note that there are four 
expressions possible for the last term. 

For n-qubit systems, we shall be particularly interested in fully separable systems which have no quantum correla- 
tions between any pair of qubits jllj. Specifically, an n-qubit density matrix p^ is fully separable on H^ n \ the tensor 
product of n two dimensional spaces, if 

P [n] =5>(a) / o (1) (a)® ■■■®p in) (a), (4) 

a 

where p^ k '(a) are qubit density matrices. If k = j, the arguments of (pf) give 



™mJ(j\p\j)(j\p\j) > max|(w|p|u)|. (5) 

J \ u 

We apply this result to the n-qubit Werner state fu^ . Define the generalized GHZ states indexed by j — Qj2- jn as 

\*±(j)) = (\j)±\j))/V2, (6) 

and the generalized Werner states as 

W ± W(s,j) = ^I n + sp ± (j) (7) 



where p ± (.7) = |^ r± (j))( 1 ^ ± (j)| and I n is the identity matrix on = <^ n Ti.2- Then choosing j judiciously in (||), 
it is easy to see that a necessary condition for the Werner state to be fully separable is that s < 1/(2™ _1 + 1). We 
shall show below that this condition is also sufficient by using the spin basis representation to give a fully separable 
expression of these states when s < l/(2 n_1 + 1) p| , ^2[ . 

We can also apply (||) to a convex set of n-qubit states which includes the generalized Werner states. Let Z?M 
denote the set of density matrices on Ti^of the form 

Jm 3m 

P (t) = Y^(tf P + (j)+tjp-(j)) £>t + t j) = i (8) 

3=0 3=0 



where if > and j m = 01 ... 1. The only non-zero elements of these density matrices are on the main positive and 



negative diagonals when the matrix is expressed in the computational basis (^) n {|0), |1}}, thus the Peres condition 
is equivalent to (||). In addition to the generalized Werner densities, 2?'"' also contains the states invariant with 
respect to depolarization: t + (j) = t~(j) for all j ^ 0. Recall that depolarization is carried out by averaging over the 
application of an arbitrary rotation exp{i<j) r a z ) to each qubit, subject to ^2 r cj> r — 2ir, followed by spin-flip of all the 
qubits This subset is the set of density matrices with non-zero entries appearing only on the main diagonal 

and in the upper and lower corners. Note that the set {/^(j)} is the set of extreme points of 2?'™', and for states in 
X)M. If p[ n ' denotes the convex subset of fully separable states in T>^ n \ a necessary condition for p(t) S is 



min 

3 



dn(t++tj)>m a x\t+-t-\. (9) 



We proved the necessary condition (^) in the computational basis; however, to find sufficient conditions we shall 
work in the spin basis. To see why this is appropriate, we first prove a useful sufficiency condition which is expressed 
entirely in terms of the Pauli matrices. 

Theorem 1 Let M n be a set of n unit vectors, M n = {mi, • • • , m„} and define the usual scalar product of the Pauli 
matrices with m 7 tr m = a ■ m. Then the density matrix 

p ± {M n ) = -^(oo ® • •• ® «r ± o- mi <g> • •• <g> CT m J (10) 

on T-v- n ' is fully separable. 
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This is easily proved using induction. The densities P ± (m) = (cro ± c m )/2 on TL2 are projections and are trivially 
separable. Suppose that p ± (M„_i) is fully separable on 7v n ^ 1 '. Then 



cr o 0[P + (m„) + P (m„)]±(T mi 



o- mn _j [P + (m„) - P (m„)]) 



= -^(M^i) ® P+(m„) + p T (M„_!) ® P-(m„)], 



(11) 



where all the upper signs and all the lower signs go together, is completely separable on TH\ n \ completing the proof. 
This particular set of separable states has the property that the only non-zero (classical) correlation is among n Pauli 
spin matrices. For n = 2 this corresponds to the case studied in reference p5|, p — (cro ® °o + Tjk<Jj ® Ofc)/4, and 

j,k 

in our case the T matrix is of rank one and T = mi ® m.2. 

To express a density in the spin basis, given its definition in the computational basis, we need a change of basis 
formula to relate the coefficients in the two bases. This is a standard exercise, but we do it in a non-conventional way 
to emphasize its structure as a Fourier transform. On the Hilbcrt space H.2 introduce the operators E a ,b — 
where a, b = 0, 1, and index the Pauli matrices using binary notation, so I2 = cr 00 , a x = croii e ^c. Both of these sets 
form a basis of operators on the space of qubit states. Write the spin basis S and the adjusted basis A in 2 x 2 arrays 



(S) 



coo 001 

(Til IfTio 



(A) 



H 



1 1 

1 -1 



(12) 



where H is the Hadamard matrix and the elements of (S) are the (real) Pauli matrices. Note that Aj t k — Pj,j©fc and 
Sj,k = Cj,j®k- In addition to the unconventional labelling in (|l2|) , it is necessary to work with real Pauli matrices, 
which is why the factor of i appears. It is then easy to check that Sj t k = J2r H 0> r ) A r ,k, which we record as 



(S) = H ■ (A). 



(13) 



The sets {Sj^/V^} and {Aj^} each form an orthonormal basis of operators on H2 where we use the trace inner 
product (B,C) — tr{B^C). Therefore, for any p 



j,k j,k 

where Sj^k = t r {Sj k p), aj t k = tr(Aj k p) = Pj,k®j and it follows that 

(s) = H.(a). 



(14) 



(15) 



The use of the Hadamard matrix allows an easy generalization to tensor product spaces. Define A^j, and S^l in 
the usual way: S^l ~ Sj lt k 1 <%> • • • ® Sj n ,k„ so, for example, = /„. It then follows easily from (S) = H ■ (A) that 

{S [n] ) = hW{AW). (16) 

The two sets of 2 2 ™ operators {S^l/V2™} and each form an orthonormal basis on the Hilbert space 

£(P>1) of 

linear operators acting on TC^ where the inner product on C is (B, C) = tr{B^C). Therefore, we can express 

an arbitrary n-qubit density matrix in the form 



'j.k 



43 



(17) 



3,k 



with [s] = H^[a], the analogue of (|15|) . As an example of the notation, both the matrices [s] and [al for a density 
matrix in the set defined in (pj) have zeros everwhere except in the first and last columns. Equation (U3) can also be 
interpreted as a two dimensional Fourier transform of the computational basis which defines the spin basis, and there 
is a natural generalization to d-dimensions using finite Fourier transforms Jig ], 

As a first application we use the spin representation to show the density matrices (0) are fully separable for 
s = 1/(2" _1 + 1). Consider VF + ^ (s, 0), but the result is independent of which j-state we choose. In terms of the 
adjusted basis, 
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The first two terms in the brackets are diagonal projections and are therefore fully separable. We write the last two 
terms in the spin coordinates. The only non-zero spin coefficients are in the last column, and 

.,,, = |(l + (-!)«»), 

where we have used ffj^ = (— \y Qk with j k denoting the binary scalar product. Define the set of 2 n_1 ele- 
ments Ind — {j : j 1 = ^2 r jr = Omod 2} . It follows from some easy algebra, that adding and subtracting a term 

In] 

proportional to the identity I n = S~ - gives 

Notice that j € Ind means that there are an even number of factors of S\ i = ia„ in S.l. so that S^ n ~ is Hermitian. 
The reason for adding and subtracting the identity is that (M) shows each term in the summation on the right is fully 
separable. To guarantee that W + ^(s,0) is a density matrix, the coefficient of the first term must be non-negative, 
forcing s < I/(2™ _1 + I) and concluding the proof that W^ n Hs,j) is fully separable if and only if s < I/(2™ _1 + I). 
This result may be compared with those obtained earlier in jllj and |Q . 

We next use the spin representation to establish a new and general sufficient condition for full separability. We 
introduce a norm on densities which is expressed in terms of the spin coefficients. 

Theorem 2 If the spin coefficients Sj^ of a density p on TC^ satisfy \\p\\i = Y]u k)^(o o) \ s j,k\ < 1> then p is fully 
separable. 

Since {—i) 3Qk is Hermitian, i 3&k Sj^ must be real. Now use (|l7|) to write 

p = a - iHix) rAl + E,, fc) , (0;0) m h (41 + ^ JQk 42) . 

where Vj t k is the sign of i 30k Sj,k- Again (|o|) applies and gives full separability for p. This guarantees that there is a 
neighborhood of the completely random state S^~/2 n in which all the densities are separable, and in particular that 
every density with \sj t k\ < 1/ (2 2 ™ — l) is fully separable, giving the analogous result in as a corollary. 

Kp = W + W(s(n),j), with s(n) = l/(2 n_1 + l), then y^ = (2™ - I)/(2™- 1 + 1). Thus condition y^ < I is sharp 
for n — 2 but may be too restrictive for larger n. One can take advantage of the special structure of a class of densities 
to obtain more refined conditions. For example (^) is also sufficient for the states inX*^ invariant with respect to 

depolarization. Consider also the following subset of DW. Let t^^ = (1 — s)/2 n ~ 1 + su^ so that ^ (tit + uj) = 1 

(recall j m = 01 . . . 1), and let /x(s) = (1 — s)S~ ra ~/2™ + sp(u). Using the same approach that was used with the Werner 

( 

densities, we find that /i(s) is fully separable provided s < 1 + 2™ 1 ^ |it+ — u . |) 

V J=6 

As our final result, we show that for any n > 2 and e > there exists a density p on Wl"! which is not fully separable 
but which has Wp^ < 1 + e. Thus, the bound of the theorem is not only the best possible in general but also the 
best possible for each value of n. As part of the proof we give necessary and sufficient conditions for full separability 
for a class of densities each of which satisfies the Peres condition. Define first the subset of with all 

diagonal elements equal: p^{t)jj — (t + (j) + t~(j))/2 — 1/2™. In the computational coordinates p^(t) is constant 
down the main diagonal and the only non-zero entries are on the main negative diagonal. Each such density matrix 
satisfies (||) and In spin coordinates has the form 

f^-^Q + ±s^. (18) 

3=6 
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If is fully separable, it can be expressed in the form (0) with p(a,jk) — (o"o + a ■ m(a, ife))/2 where m(a,jk) is a 
unit vector in the x-y plane. Let 9(a,jk) be the angle m(a,jfc) makes with the x-axis. Then from (^) the non-zero 
off-diagonal elements are 



n—x 

a r—0 

All the matrix elements of the matrices in are real, so the densities in Dc satisfy 

71—1 

T£ p (a)-cosC£(-l)H(a,j r )) = P% 

a r=0 

It is immediate from ( p(i| ) that if p?-. = 1/2™ for some j, then for all a 

n-i 

^r(-l) jr 0(a,j r ) = mod 2tt. 



(19) 



(20) 



(21) 



r=0 



We use (|2l]) to define necessary and sufficient conditions for full separability for a subset T>c of T>c ■ Included in 
T>^ are non-separable density matrices with HpI^ arbitrarily close to 1, confirming the assertion that \\p\\i < 1 cannot 
be improved for n x n densities. To illustrate the ideas with minimal notational clutter, we work with the case n = 3. 
Let pooo,m = pooi,no = 1/8- Then ( |2l| ) implies that for all a, 

9(a,l) + 6(a,2) +6(a,3) = modl-ix and 6 (a, 1) + 9 (a, 2) - 6 (a, 3) = mod 2%, 

so that 9 (a, 3) = and 9 (a, 2) = —6 (a, 1). But then it follows that a necessary condition for full separability is 

Poio,ioi = Pon, ioo = \ ^2p{a)cos(29 (a, 1)). 



Define as the set of states with the additional restrictions: if c and d satisfy —1/8 < c, d < l/8,thcnf ± (010) = |±c 
and t^ 1 (011) = h ± d. Thus the states in have the form 



1/8 

1/8 

1/8 

1/8 





0c 

1/8 

1/8 






d 

1/g 
d 





1/8 
1/8 






1/8 
1/8 
1/8 



(22) 



It follows that p (t(c, d)) is not fully separable if c ^ d. Using spin coordinates, it is easy to establish that \\p (t(c, d))\\ = 
1 + 4 |c — d\ , and thus p (t(c, d)) is fully separable if c = d. Since \\p (t(c, d))\\ can be made arbitrarily close to 1, we 
have shown \\p\\i < 1 is sharp for n = 3. The argument for larger n is similar, and we omit the details. 

Proposition 3 Let iffl denote the subset of densities in T>\^ which, in addition to being constant on the main 
diagonal, have t + — l/(2™ -1 ) and t~ = for the first (and last) 2 n ~ 2 positions on the main negative diagonal, 
= ± c on the next 2™~ 3 positions and ^ = A ± rf on the remaining positions, where —1/2" < c,d < l/2 n . 

Then every density in T> c satisfies the Peres condition and is fully separable if and only if c = d. Given e > 0, there 
exist densities in T>^ which are not fully separable and \\p\\i < 1 + e. 
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